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Abstract 

We present the mechanism of interaction of Wnt network module, which is responsible for periodic 
sometogenesis, with p53 regulatory network, which is one of the main regulators of various cellular 
functions, and switching of various oscillating states by investigating p53—Wnt model. The variation 
in Nutlin concentration in p53 regulating network drives the Wnt network module to different states, 
stabilized, damped and sustain oscillation states, and even to cycle arrest. Similarly, the change in Axin 
concentration in Wnt could able to modulate the p53 dynamics at these states. We then solve the set 
of coupled ordinary differential equations of the model using quasi steady state approximation. We, 
further, demonstrate the change of p53 and Gsk3 interaction rate, due to hypothetcal catalytic reaction 
or external stimuli, can able to regulate the dynamics of the two network modules, and even can control 
their dynamics to protect the system from cycle arrest (apoptosis). 
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Introduction 

p53, one of the largest hub in cellular network mm, is considered to be one of the most important 
key regulators of cellular metabolic pathways Oil], regulates a number of cellular functions [5] and its 
dynamics control even the fate of the cell when stress is given to the cell [^. Further, p53 suppression is 
observed in various types of cancer either due to mutations or by ambiguous expression of control systems 
like MDM2 (or HDM2 its human equivalent) [7]. It is negatively regulated by the MDM2 [5], and 
functions as an E?> ligase which mediates the proteosomic degradation of p53 mni, facilitates the diffusion 
of p53 in the nucleus m, which in turn lowers the level of p53 within the cell [T2]. In stressed cells, the 
p53 triggers the cell to cell cycle arrest forcing it to choose its fate, either to repair or apoptosis [T3lll5) . 

p53 dynamics is regulated by various signaling molecules as evident from various experimental and 
theoretical reports [16]. On of the most important inhibitor oip53 — MDM2 is nutilin m- It is a selective 
small molecule which directly binds to binding pocket of MDM2 [18], to activate p53 pathway [TSlIlQj . 
and inhibiting p53 binding to MDM2 |20j . Nutilin — 3a and 36 are the two most active enantiomorphs 
that are found to up-regulate the p53 in p53 dependent manner [21]. Further, it is also demonstrated 
that nutlin-Z induces anti-angiogenic activities in endothelial cells probably via three mechanisms; first 
by inhibiting endothelial cell migration; second by inducing cell cycle arrest; and third by increasing 
apoptotic tendency in endothelial cells [22] . Further, it was also shown that nutlin-S treatment in these 
cells leads to accumulation of p53 [23] . indicating the important impacts of nutlin-3 which interferences 
the pb‘i-MDM2 regulatory mechanism. As a consequence, it was shown that Nutilin-3 stabilizes the 
p53 dynamics and causes the activation of survival pathways |23| . 

Somitogenesis in vertebrates is periodic formation of somites (vertabrae precursors) in the anterior 
presomitic mesoderm tissue (PSM) controlled by complex gene network known as segmentation clock [24], 
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where, Notch, Wnt pathways and fibroblast growth factor (FGF) are the main components |27) . 

Further, it has been shown that these three main pathways generate rhythms of specific oscillation rela¬ 
tionships [37] , and cross-talk among them [35] as a basis of spatio-temporal self-organization of patterns 
during sometogenesis [33]. This physiological oscillations is responsible for periodic spacing of somites, 
dynamic structures and regular segmented development in vertebrates |29| . In some studies in normal 
(stress) and cancerous cells, it has been reported that p53 and Wnt modules are being coupled and 
cross-talk between them via different intermediate proteins or genes or signaling molecules [301132] . The 
studies also reported that p53 network suppresses the Wnt signaling cascade [31]. The Wnt signaling 
plays a pivotal role in determination of cell fate, which may probably through p53 interaction Eiiisa, 
but still it is not clear yet. Further, several other works reported that the intercellular signaling of Wnt 
depends on the /3-catenin and glycogen synthase kinase-3 {GSK-3) [33]. Moreover, p53 can regulate the 
dynamics of Axin2 and GSK-3 [33]. Recently, it is also reported that the miR — 34 family, which is 
directly transcribed by the p53, links p53 with Wnt signaling [33], in this process, and a set of GSK-3- 
related Wnt genes, are directly targeted by p53 and mii?-34. These observations clearly suggest the close 
connection between the p53 functionalities and Wnt in stress and cancer cells [30l|33]. However, the way 
how p53 interacts with Wnt during sometogenesis and its impact on somite organization/reogranization 
at molecular level is still an open question. 

GSK-3, which is an important coupling molecule between p53 and Wnt pathways, is widely expressed 
in broad range of cellular processes and being a Ser/Thr kinase it is a multifunctional protein |351l36j . It 
is active under resting condition and act as a core regulator in various disease pathways like cancer |35) . 
Being a kinase, it confers selectivity and substrate specificity [35]. It also largely acts as a phosphorylating 
and an important signaling agent for the degradation of large number of proteins, such as /3-catenin in 
Wnt pathway m- Further, it is a central player in Wnt signaling pathway recruited with the Axin2 at 
the Frizzle receptor for the activation of Wnt signaling [38]. In absence of Wnt signals GSK-3 shows its 
activity in the protein destruction complex with Axin2, APC (adenomatous polyposis coli) and other 
partners mediating the destruction and phosphorylation of /3-catenin [39]. GSK-3 also mediates p53 
dependent apoptosis signaling by suppressing p53-AKT pathway. The role GSK — 3 as a mediator of p53 
and Wnt cross-talk is not fully understood. The impact of other Wnt regulators in Wnt pathway and 
how they regulate p53 dynamics, and vice versa, are needed to be investigated in order to understand 
how segmentation clock works. 

We study the coupling of Wnt oscillator which is one of the three important clocks in sometogenesis in 
vertabrates and p53 regulatory network to understand how does these two oscillators regulates each other 
and process information during somite formation and pattern organization. Since p53 interferes various 
other cellular functions, it will be interesting to understand its regulating impacts on segmentation clock 
and vice versa. Presently we study these phenomena by modeling this p53-Wnt two oscillator model 
with updated interaction between them from various experimental reports and how do they regulate each 
other. We organize our work as follows. The detailed description of the model and techniques is provided 
in materials and methods section. It is followed by mathematical treatment of the model, the numerical 
simulation results and its discussion in the next section. We then draw come conclusions based on the 
results we obtained and discussion on it. 


Materials and Methods 

A. Coupling p53-Wnt Oscillators Model 

The Wnt signaling pathway is authorized by /3-catenin regulation via various intermediate interaction 
steps (Fig. 1) [301141] . In this signal processing, Wnt promotes the synthesis of Dishevelled (Dsh) protein, 
and then inhibits GSK3 protein kinase m- GSK3, then interacts with Axin (with rates fci and k 2 for 
forward and backward reactions as in Fig.l) in the presence of /3-catenin to form transient complex which 
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Figure 1. The schematic diagram of biochemical network of p53 and Wnt modules and their cross talk 
via GSK3. 


in turn phosphorylates ^-catenin into /3-catenin-P (with rate rate ko) [27l|4T]. The /3-catenin-P degraded 
quickly with rate ^ 14 , on the other hand it is dephosphorylated with rate to form freed /3-catenin and 
then degraded with a rate ^43 [32]. The free /3-catenin is then transported to nucleus (/3-catenin-n) with 
rate ^ 4 , on the other hand /3-catenin-n is transpoted from nucleus to cytoplasm with another rate k^. 
Then /3-catenin-n promotes the transcription of set of target genes including mAxin (mRNA of Axin, 
specifically homolog Axin2 in PSM tissue) at the rate fcis, with rate kie, and mAxin also degrades with a 
rate kis [43] . On the other hand mAxin translates to Axin with a rate Axin protein which also degrades 
with a rate /c 23 - This mechanism of Axin is just like a negative feedback loop which drives the oscillations 
in the network components. Thus Axin2 is the essential target of the Wnt pathway because it inhibits 
the Wnt signals by degrading the /3-catenin by forming a negative feedback loop [44] . 

The Wnt pathway communicates with the p53 Oscillator [45ti47] through GSK3 by formation of a 
binary complex p53 — GSK3 with the free available p53 at the rate k^g. Then this complex p53 — GSK3 
degrades at the rate ^49 to the free p53 and GSK3 proteins [ISj. p53 and MDM2 regulate each other 
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with feedback mechanism where formation and decay of various intermediate complexes are involved with 
different rates (Fig. 1) as described in Table 1 and 2 [3^. The freely synthesized Nutilin directly interacts 
with MDM2 with a rate k^Q to form a complex NutMDM2, and then this complex dissociates to release 
back free Nutlin and MDM2. Then Nutilin indirectly interacts with p53 via MDM2 by affecting the 
feedback mechanism between p53 and MDM2. The Nutilin then degrades with rate constant fcag. Thus 
the activated p53 can now go to the nucleus and initiate the transcription of a number of target genes 
including MDM2. MDM2 is essential in controlling the concentration of p53 through a positive feedback 
loop mechanism by making complex |491I50) . 


Mathematical model of the network: Quasi steady state approximation 

The stress p53 — Wnt regulating network we study is defined by = 13 molecular species (Table 1) 
corresponding to the reaction network description provided in Table 2. The state of the system at any 
instant of time t is given by the state vector, x{t) = {xi,X 2 , ■ ■ ■ , where, iV = 13 and T is the 

transpose of the vector. By considering feedback mechanism of in p53 and Wnt oscillators and coupling 
reaction channels of the two oscillators, we could able to reach the following coupled ordinary differential 
equations (ODE), 
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The coupled ODEs (equation ([T|)) are very difficult to solve analytically. However, one can get ap¬ 
proximate analytical solution of these equations by using quasi-steady state approximation [5TJ[S2]- In 
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general, any biochemical reactions network involves two basic types of reaction, namely slow and fast 
reactions m- Therefore, the fV = 13 variables in the system can be divided into sets of slow and fast 
variables respectively. If x^{t) = {xi,X 2 , ■ ■ ■ ,xi)'^ and = {xi+i,xi+ 2 , ■ ■ ■ ,xn)'^ are slow and fast 

variable vectors, then x{t) = (^x^,x^) . Then from equation ((I|)-(14) along with Table I, we have, 


F(t) 
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(15) 


Since the rate of complex formation is fast, and after this fast coirrplex formation, the fast variables 
immediately retain steady state (equilibrium). Then using Henri-Michaelis-Menten-Briggs-Haldane ap¬ 
proximation |53j . one can take quite fair assumption that the ODEs of variables of complex molecular 
species reach steady state equilibrium quite fast compared with the time evolution of slow variables |51) . 
Then one can straight forward put, w 0 [511113] • Now, the number of coupled ODE (|T]) reduces to 
the following, 
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and the x-^ become the following steady state values. 
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(17) 


The thirteen ordinary differential equations are now reduced to six ordinary differential equations. This 
allows to simplify the complex mathematical model to get approximate solutions of the system numerically 
saving computational cost or analytically fron the reduced system if possible. 

We used standard Runge-Kutta method (order 4) of numerical integration to simulate equation (re- 
fode) to find the solution of the variables listed in Table 1 for the parameter values given in Table 2. We 
then analysed the constructed mathematical model to get possible approximate analytical solutions of 
the variables (slow variables) using quasi-steady state approximation. 


Results and Discussion 

The cross talk between p53 and Wnt via Gsk3 subjected to different stress conditions induced by nutlin, 
as well as Axin concentrations available in the system. In order to understand how they regulate each 
other, we numerically simulated the coupled ordinary differential equations © with the parameters listed 
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Figure 2. The dynamics of Axin and /3-catenin for different values of ^35 i.e. 0.06, 0.18, 0.24, 0.6 and 
1.2 for fixed vale of ^22 = 0.02. The right hand panels are the two dimensional plots the Axin and 
/3-catenin with nutlin for the same set of parameter values showing different state behaviors. 


in Table 2. The simulation results for the variable in the model are presented and compared to understand 
the switching of different oscillating states induced by different stress conditions. The reason is that these 
states of the oscillations in Wnt network are responsible for different behavior of regular somite formation, 
and many other periodic phenomena during sometogenesis. The roles of Axin2, Gsk3, , 8 -catenin and 
Nutilin on the p53 network dynamics is discussed in the context of the model we presented. Further, 
the signal processing between the two coupled network modules via the Gsk3 during sometogenesis is 
studied under the stress p53 generating effect due to Nutilin-Sa. On the other hand, we also studied the 
regulation of p53 by Wnt due to available Axin concentration in the system. 

Driving Wnt oscillating states by p53 

The p53 in p53 regulatory network can be induced stress by the available concentration of stress inducer 
molecular species, Nutlin —Sa. The concentration level of Nutlin —3a is proportional to the creation rate 
constant of Nutlin — 3a, /cas (Fig. 1). Therefore, computationally we can able to monitor the amount of 
stress induced to p53 by changing the value of parameter k^^ and supervising the dynamical behavior of 
p53 (Fig. 2 and 3). 

The dynamics of p53 (Fig. 2 upper left panels) for small values of ^35 (fc35(0.07) show single spike due 
to sudden stress, and due to small ^35 values, the dynamics become stabilized indicating normal nature 
of the system. Then, as the value of k^^ increases, p53 dynamics show damped oscillation indicating 
the inducing oscillation for certain interval of time {t ss [0,4000] minutes for ^35 = 0.18), and then its 
dynamics become stabilized. The time interval of damped oscillation regime depends on the value of ^35 
and increases as ^ 35 . This indicates that if the stress given to the system is small, the system first will go 
to the stress or excited state (indicated by oscillatory behavior), repair back the changes in the system 
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Figure 3. The similar plots for p53 and GSK3 for the same set of values as in Figure 2. 


and come back to the normal condition. For sufficiently large values of (0.23 < ^35 < 0.55) the p53 
dynamics become sustain oscillation, where the responsed stress is optimal and sustains the stress for 
t —> 00. Further increase in the value of /C 35 (^35 > 0.55) force the sustain oscillation of p53 to damped 
oscillation whose time interval of damped oscillation decreases as ^35 increases. This indicates that excess 
stress in p53 due to Nutlin may become toxic to the system reflected in p53 dynamics. If we increase the 
value of ^35 (^35 > 1.3) the p53 dynamics become stabilized. This may probably the case of apoptosis 
of the system. The results obtained are in consistent with the experimental observations which indicates 
that acetylation of p53 is responsible for its activation and stabilization [54] . If we further increase the 
value offc 35 ,p 53 activation decreases maintainingp53 stability, but at higher values (<35). 

The two dimensional plots between p53 and Nutlin (Fig. 2 upper right panels) for different values 
of fc 35 show different behaviors of the system, namely stable fixed point, damped oscillation towards a 
fixed point, sustain oscillation, then damped and stabilized fixed point. This behaviors show the same 
behaviors as shown in p53 dynamics for different values of 

Similarly, we study the dynamics of Gsk3 as a function of ^35 (Fig. 2 lower left panels) for corre¬ 
sponding values which we took for the case of p53. In this case also we found the similar behaviors as 
we got in p53 case (Fig. 2), namely, stabilized, damped, sustain, then damped and stabilized states for 
the corresponding values of ^35 taken in p53 case. The two dimensional plots of GsfcS and Nutlin (Fig. 
2 lower right hand panels) show the evidence of the existence of the mentioned states. 

Now we study the dynamics of Axin (Fig. 3 upper left panel) imparted by stress p53 via Gsk3 induced 
by Nutlin concentration in the system. The signal of the stress p53 is processed by Axin via Gsfc3 and 
reflected the same states in Axin dynamics as obtained in p53. The Axin and /3 — Gatenin dynamics 
as a function of ^ 35 , as well as two dimensional plots of these two molecules with Nutlin are similar 
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Figure 4. Plots of p53, Axin and /3-catenin as a function of time for ^22 values 0.0002, 0.0003, 0.01, 
0.025 and 0.08, for fixed values of ^ 35 . The upper right panels shows the two dimensional plots of p53 
and /3-catenin with Axin for same set of ^22 and ^35 parameter values. 


qualitatively with the dynamics of stress p53 and GsfcS (Fig. 2 and 3). This means that depending on 
the values of ^35 the stress p53 can able to generate and arrest the Axin and j3 — Catenin cycles. Thus the 
dynamics of Axin and /3 — Catenin can be regulated by p53 and probably can control the Wnt oscillator. 
This means that p53 regulation has strong impact on Wnt oscillation, and regulate various mechanisms 
during sometogenesis. 

Dynamics of p53 regulated by Wnt 

The activation of Axin occurs after the Gsk initiation from the signal received at the membrane. The 
Nutilin synthesis rate (^ 35 ) is kept constant in this case (^35 = 0.36), giving a slight stress to the system. 
In order to understand how does Wnt oscillator regulate p53, we allow to change Axin synthesis rate ^22 
and see the dynamical behavior of p53 induced by Axin, which is one of the most important molecular 
species in Wnt network module. The dynamics of the Axin and p53 for different values of k 22 (Fig. 4). 
For small values of ^22 (^22 < 0.0002), Axin dynamics maintains lowest constant level, showing stabilized 
state of the Axin dynamics (Fig. 4 lower left panels). The increase in the value of ^22 allows the Axin 
dynamics to switch from stabilized state to damped state for certain range of time and then become 
stabilized. The co-existence of these two states (damped and stabilized states) takes place for a certain 
range on ^22 (0.0003 < ^22 < 0.008). If we increase the value of ^ 22 , the duration of damped oscillation 
also increases. The switching of damped state to sustain oscillation state takes place for certain range of 
^22 , 0.00 87 < ^22 < 0.015. Further increase in ^22 force the sustain oscillation state to damped oscillation 
state again, where the time interval of damped oscillation becomes smaller as ^22 increases. If we increase 
^22 value, then Axin dynamics landed to stabilized state. 
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Figure 5. The amplitudes of p53 (x^), /3-catenin {xf), Gsk3 (xg ) and Mdm2 (xg ) as a function of ^35 
for various fcgg values 0.030, 0.035, 0.040, 0.045 and 0.050 (upper four panels). The lower four panels are 
the maxima of p53 (xf^), /3-catenin (xf^), Gsk3 (x|^) and p53 (x^) for the set of parameter values. 


We then show the dynamics of /3— Catenin as a function of /C 22 (Fig. 4) which indicates the switching 
of the three different states driven by fc 22 ■ In this case we could able to find sustain oscillation very quickly 
indicating that the ,5—Catenin can be quickly switch to stress condition and easily. Similarly, /3—Catenin 
can reach second stabilized state (probably cycle arrest) quickly and easily as compared to other molecular 
species. The dynamical behavior of ,5—Catenin is supported by two dimensional plots of /?—Catenin and 
Axin (Fig. 3 middle right panels) which exhibit fixed point oscillations (corresponding to two stabilized 
states), attractor towards the stable fixed points (corresponding to two damped oscillations), limit cycle 
(corresponding to sustain oscillation). 

Now the dynamics of p53 as a function of ^22 (Fig. 4 upper left panel) show that the patterns in 
the dynamics of molecular species in Wnt oscillator are well captured in the p53 dynamics. The p53 
dynamics exhibit the corresponding states, namely, stabilized, damped and sustain oscillation states, as 
observed in Axin and /3—Catenin as a function of ^22 with similar patterns and dynamics. 

The important dynamical behaviors those have been captured in p53 dynamics due to change of 
states in molecular species in Wnt oscillator are synchronous in states and time periods of dynamics. If 
we look at the dynamics of Axin and /3—Catenin as a function of /c 22 j "we can easily observe two important 
changes, ( 1 ) change in time period of oscillations of the molecular species, and ( 2 ) change in the pattern 
of states. The increase in ^22 force to decrease the time period of the Axin and /?—Catenin dynamics. At 
the same time the switching of the various states takes place. These changes in the Axin and /?—Catenin 
are systematically and synchronously well captured in qualitative sense in the p53 dynamics (Fig. 4) 
during the communication of the two p53 and Wnt oscillators via Gski. 
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Figure 6. Similar plots of p53 (xi), /3-catenin (is), Gsk3 (xe) and p53 (xr) as a function of k 22 for the 
same ^39 values as in Figure 5. 


Amplitude death driven by Gsk3 and p53 interaction 

We now study the impact of interaction of GskS and p53 (^ 49 ) on Wnt and p53 cross talk driven by 
Nutlin concentration level in the system (Fig. 5). The rate of interaction of p53 and Gsk3 can modulated 
by external catalytic action or external stimuli regulating the interaction in the system. 

Wnt cycle arrest induced by Nutlin 

We calculated amplitudes of four molecular species Axin, /?—Catenin, Gsk3 and Mdm2 as a function 
of ^35 for five different values of ^46 (Fig. 5 upper two rows). The amplitudes of the molecular species 
for every values of ^42 are calculated by averaging the amplitudes within the interval of time [100,1500] 
minutes. From the results we observe three different regimes for each value of ^ 45 : first amplitude death 
(AD) regime (for small ^ 35 ); second finite amplitude or oscillation regime (OR) (for moderate ^ 35 ); and 
third amplitude death regime again (for large fc 35 ). The first AD regime may correspond to normal state 
(stable state) of the system, and second AD regime may correspond to apoptosis (cycle arrest) of the 
system. Further, we calculated maximum values of molecular species Axin, /3—Catenin, GsfcS and p53 
for the correponding values taken above (Fig. 6 lower two rows of plots). The patterns obtained in these 
plots are similar to the patterns of the amplitude plots of the corresponding molecular species. 

In OR, the amplitude is finite and increases monotonically as a function of ^ 42 , then reaches a 
maximum amplitude, after which amplitude decreases slowly and then decreases monotonically. This is 
the regime where the system is in activated or stress condition. It is also to be noted that as fe 46 increases 
the amplitude also increases and the range of activated regime also increases. This means that increasing 
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^46 force the system to get activated quickly and protects the system from next stabilized regime (cycle 
arrest or apoptosis regime). Therefore, even though increase in Nutlin concentration drives the system to 
apoptosis regime quickly, one can modulate the interaction of p53 and Gsk3 to protect the system from 
apoptosis upto significant range of Nutlin concentration level, so that the system can repair back from 
toxic due to excess Nutlin to come back to normal condition. 

p53 cycle arrest induced by Axin 

Now we study the impact of p53 and GskZ interaction on the condition of activation of molecular species 
of both the p53 and Wnt oscillators. We calculated the amplitudes of these molecular species as a function 
of Axin concentration level in the system (^ 22 ) (Fig. 6 first two rows) for the same parameter values as 
we had taken in the case of calculation of amplitudes of these molecular species with the variation of 
Nutlin concentration in the system. For the same parameter values we also calculated the maxima of 
these molecular species (Fig. 6 lower two rows). 

The change in ^22 drives the amplitudes of all the molecular species to different oscillation states 
(AD, OR, AD) as we have obtain in the case of Wnt cycle arrest by ^42 variation. The increase in 
Axin concentration level in the system forces amplitudes of the molecular species variables to pass the 
trajectories of first AD (normal condition), then OR (stress condition) and then second AD (cycle arrest 
condition). However, in this case the amplitudes of the variables do not change much as a function fc 22 , 
but significant change in the range of second AD as a function of k^Q. The results indicate the Axin 
induced p53 cycle arrest, however, interaction of p53 and Gsk3 can able to protect the system from 
apoptosis upto significant range of ^22 so that the system can repair back to its normal condition. 

Steady state solution of the system 

Since the direct and exact analytical solution of the set of coupled nonlinear differential equations o, 
we used quasi steady state approximation [5T1[52] to obtain approximation solution of the system. Using 
steady state values of the fast variables using dm, putting back to the ODE of xi{t) in equation (|T 6 )l 
and taking ^ (^ ~ ■^) ^ ^^i) we get the following approximate solutions of 

Xi, 


Xi{t) « —+Cie-“=‘ 

U2 

where, Ci is a constant, and ui, U 2 are also given by, 


(18) 


Ml 

U2 


k26 


(^22X2 + ^ 2^8 


142^2 

klk2e 



^1^26^40 2;i3 
(^41X33 — fc 29 X 9 ) 


(19) 


The equation (IT51) shows that Xi becomes constant (steady state) as t —>■ 00 . However, xi depends on 
various other steady state variables as shown in equation (1181) and m- 

We then solve ODE of X 3 (t) in equation (ITCl) for X 3 , and the result is given by, 


^3(t) « ^ 

M3 


( 20 ) 


where, G 2 is a constant. Here, in this solution, we can get that X 3 stabilizes to a constant value at 
sufficiently large time limit {t —^ 00 ). 
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Now, the approximate solution of Xq can be obtained by solving ODE of xq in equation dT51) . using 
equation (fT^ and (IT71) . as given below, 


x^it) ^ exp {at — be 



^a/u2p 

U 2 J 



( 21 ) 


where, C 3 is a constant. Similarly, U 3 = ^ 2 ^ 8^261 a = and b = gj^g^ constants. 

The molecular species Xn can be solved directly from equation (1161) and using equation (1171) . The 
result is given by, 

xii(t) + ^ (22) 

kss 

where, C 4 is a constant. Further, using the differential equation in 0:7 in equation dSD and equation (1221) 
along with equation (II3, we solve for 0:7 given by. 


X 7 {t) 


U4 r 
kss J 


y-Ci 


dy + Cs 


Ca + 1 ^ 6^384 


+C 5 

Ca + ff 6^=38* 


(23) 


where, ua = ^31 + ^ 34 x^ 0 , U 5 = (^ 33 ^ 37 * 12 ) and C 5 is a constant. Similarly, the variable Xg can be 

solved using equations dSD, (EH) and (II3 as follows. 


Xs{t) 


Cfi- 


k27xl 

^38 


Jy 


- i+i 




dy 


-UGt-\-^e-^38^ 

e «38 




(s) 


e"“'^‘+'=38 




(24) 


where, Uq = M 7 =^ 36 C'i and Cq are constants. E is the error function. 


Conclusion 

The cross talk between p53 and Wnt oscillators can influence each other such that the dynamics one 
oscillator can regulate the dynamics of the other and vice versa. The concentration level of Nutlin can 
drive the p53 network module to different oscillating states corresponding to different stress states of the 
system. This changes in the oscillating states can be well communicated to Wnt network module via 
Gsk3 and could able to regulate the Wnt functionalities. Further, this changes can well interfere and 
give impact on the process of sometogenesis of Wnt network module. On the other hand, the available 
concentration level of Axin in Wnt network module can also influence the dynamics of molecular species 
variables in Wnt network module at different stress states. Then the information of these changes is well 
process by the p53 network module via GsfcS and capture the changes accordingly. This means that the 
p53 and Wnt network modules can cross talk between them and regulate each other due to changes in 
each oscillator by processing information between them. However, excess changes in each network module 
can act as toxic to each module and arrest the cycle of the other, and vice versa. 

We then demonstrate that any other changes in the regulation of the network module, may probably 
by catalytic reaction or external stimuli which alters the rate of the reaction of a particular reaction, can 
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probably try to protect the system from apoptosis upto certain range of system’s condition within which 
the system can repair back the changes to come back to normal condition. We attempt to numerically 
demonstrate this phenomena by allowing the change in the rate constant of p53 and GskS interaction and 
clearly give a chance to both the p53 and Wnt network modules upto certain range of excess concentration 
levels of Nutlin and Axin to protect themselves from cycle arrest, namely apoptosis. This means that one 
important role of catalytic reaction and external stimuli (which do not affect the reactions and network 
but rate of reactions) is to protect the system from apoptosis or system’s failure by giving a chance of 
repair the failure or defect by itself within the system automatically. Therefore, we need to find out 
various such reactions and stimuli in cellular network in order to control and prevent the cell from cell 
death, due to cellular breakdown by diseases or other external mechanisms, by regulating these important 
reactions and stimuli. 
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Table 1-List of molecular species and their intial concentration 


S.No. 

Species Name 

Description 

Notation 

Intial Concen- 
tration(nM) 

1 . 

Axin2 

Intial Concentration Axin2 protein 

Xi 

10.0 

2 . 

Aa;m2 mRNA 

Intial Concentration Aa;m2 mRNA 

X2 

0.1 

3. 

/3-catenin 

UnPhosphorylated /3-catenin 

X3 

0.1 

4. 

/3-catenin p 

Phosphorylated /3-catenin 

X4 

0.001 

5. 

/3-catenin n 

Nuclear /3-catenin 

X5 

0.1 

6 . 

Gsk3 

Gsk3 protien 

xe 

5.0 

7. 

p53 

unbound p53 protein 

X7 

0.022 

8 . 

MDM2 

Unbound MDM2 protien 

Xs 

0.035 

9. 

MDM2 mRNA 

MDM2 Messsenger mRNA 

Xq 

0.01 

10 . 

p53.MDM2 

MDM2 and p53 complex 

Xio 

0.01 

11 . 

Nutlin 

Unbound Nutlin 

Xll 

0.0425 

12 . 

Nutlin.ML»M2 

Nutlin MDM2 complex 

Xl2 

0.01 

13. 

p53-Gsk3 

p53 and Gsfc3 complex 

Xl3 

0.01 
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Table 2 List of parameter 


S.No. 

Notation 

Description 

Values 

Reference 

1 . 

Ki 

Rate constant for binding of Gsk3 to Axin2 

0.2SnMmin~^ 

m 

2 . 

K 2 

Rate constant for dissociation of Gsk3-Axin2 

0Amin~^ 

m 

3. 

Ks 

Rate of /3-catenin synthesis 

0.087nAImin~ 


4. 

K 4 

Rate of /3-catenin entry into the nucleus 

0.7min~^ 

m 

5. 

K 5 

Rate of /3-catenin exit from the nucleus 

1.5min~^ 

m 

6 . 

Kq 

Rate of phosphorylation of ^-catenin by the Gsk3 

5.08nMmin~^ 

m 

7. 

Kr 

Concentration of Dishevelled(Dsh)protein 

2.0nM 

m 

8 . 

Ks 

Total Gsk3 Concentration 

3.0nM 

m 

9. 

Kg 

Rate of inhibition by Dsh for/3-catenin_p by the 
Axin2JGsk'i complex 

O.bnM 

m 

10 . 

Kio 

Rate constant for /3-catenin_p by the Axin2_Gsk3 
complex 

0.28nM 

m 

11 . 

Kii 

Rate constant of dephosphorylation of /3-catenin 

1.0nMmin~l 

m 

12 . 

Ki2 

Maximun rate constant for /3-catenin phosphora- 
tion 

Q.OOSnM 

m 

13. 

Ki3 

Rate constant for degradation of unphosphory- 
lated /3-catenin 

0 

m 

14. 

Ku 

Rate constant for degradation of phosphorylated 
/3-catenin 

7.062mm—1 

m 

15. 

Ki5 

Rate of transcription of the Axin2 gene 

0.06nMmin~^ 

m 

16. 

Kie 

Rate of transcription of the Axin2 gene induced 
by nuclear /3-catenin 

1.6AnMmin~^ 

m 

17. 

Ki7 

Rate for induction by nuclear /3-catenin of Axin2 
gene trascription 

0.7nM 

EU 

18. 

Kis 

Maximum Rate of degradation of Aa;m2 mRNA 

0.8nMmin~^ 

m 

19. 

Ki9 

Rate constant for degradation of Aa:m2 mRNA 

0A8nM 

m 

20 . 

K20 

Rate of transcrioption of Axin2 gene induced by 
transcription factor 

0.5nM'min~^ 

m 

21 . 

K21 

rate for induction by transcription factor of 
Axin2 gene transcription 

O.ObnM 

m 

22 . 

K22 

Rate of synthesis of Axin2 protein 

0.02min~^ 

m 

23. 

K23 

Maximum rate of degradtion of Axin2 protein 

0.6nMmin~'^ 

m 

24. 

K24 

Michaelis rate for degradtion of Axin2 protein 

OmnM 

m 

25. 

K25 

Hill coefficients 

2.0 

m 

26. 

^26 

Scaling factor for Wnt oscillator 

1.5 

m 

27. 

^27 

MDM2 synthesis 

0.297mm“^ 

m 

28. 

^28 

MDM2 transcripton 

0.006mm“^ 

m 

29. 

feg 

MDM 2 mRN A degradation 

0.006mm“^ 

m 

30. 

kso 

MDM2 degradation 

0.2598mm“^ 

m 

31. 

ksi 

p53 synthesis 

4.68mm“^ 

m 

32. 

^32 

p53 degradation 

0.0495mm“^ 

m 

33. 

kss 

p53-MDM2 binding 

0.693mm“^ 

m 
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Table 2 Continue... 


S.No. 

Notation 

Description 

Values 

Reference 

34. 

k34 

p53J\IDM2 dissociation 

0.00693mm“^ 

m 

35. 

k35 

Nutilin formation 

O.OOlmm”^ 

Estimated 

36. 

^36 

Nutilin MDM2 formation 

0 .012mm“^ 

Estimated 

37. 

k37 

Nutilin MDM2 degradation 

0.03771*71“^ 

Estimated 

38. 

00 

CO 

Nutilin degaradation 

0.06r?7777“^ 

Estimated 

39. 

CO 

p53-Gsk3 complex formation 

0.04777777“^ 

m 

40. 

kio 

p53-Gsk3 complex degradation 

0.12777*77“^ 

m 

41. 

kii 

MDM2_mRNA synthesis 

0.042777*77“^ 

m 
















